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ABSTRACT 



We analyze the lightcurve of the microlensing event OGLE-2003-BLG- 
175/MOA-2003-BLG-45 and show that it has two properties that, when com- 
bined with future high resolution astrometry, could lead to a direct, accurate 
measurement of the lens mass. First, the lightcurve shows clear signs of distor- 
tion due to the Earth's accelerated motion, which yields a measurement of the 
projected Einstein radius f E . Second, from precise astrometric measurements, 
we show that the blended light in the event is coincident with the microlensed 
source to within about 15 mas. This argues strongly that this blended light is 
the lens and hence opens the possibility of directly measuring the lens-source 
relative proper motion /j, rel and so the mass M = (c 2 /4G)/i rel t E f E , where t E is 
the measured Einstein timescale. While the lightcurve-based measurement of f E 
is, by itself, severely degenerate, we show that this degeneracy can be completely 
resolved by measuring the direction of proper motion /x rel . 

Subject headings: astrometry - gravitational lensing - stars: fundamental pa- 
rameters (masses) 



1. Introduction 

When microlensing experiments were initiated more than a decade ago (Alcock et al. 
1993; Aubourg et al. 1993; Udalski et al. 1993), there was no expectation that the individual 
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lens masses could be determined to much better than an order of magnitude. The only 
routinely observable parameter, the Einstein timescale t E , is related in a complicated way 
to the mass M and two other parameters, the lens-source relative parallax, 7r re i, and relative 
proper motion /i Te i. 

* = (i) 

where 

9 E = ^kMh xcX (2) 

is the angular radius of the Einstein ring and k = 4G/(c 2 AU) ~ 8.14mas/M . In principle, 
therefore, a measurement of 9 E and 7r rcl would lead to a determination of lens mass (Refsdal 
1964). However, since neither 7r rc i nor /i rcl are usually known for microlensing events, one 
can generally obtain only a rough estimate of the lens mass based on statistical inferences 
from the distance and velocity distributions of the lens and source populations. 

The motion of the Earth in its orbit produces a distortion in the observed lightcurve 
from that of the simple heliocentric case. The magnitude of this distortion is proportional 
to the size of the projected Einstein radius r E relative to the size of the Earth's orbit. This 
ratio, 7Te = AU/fg, is commonly called the microlens parallax, from the similarity in its 
definition to astrometric parallax. As shown in Gould (2000), n E 9 E = 7r rcl , and therefore 



Gould (1992) pointed out that individual lens masses could be determined provided 
that #e and tc e were simultaneously measured for the same event, and he suggested some 
methods for measuring each 1 . If successfully carried out, microlensing would join only a 
handful of other methods for directly measuring stellar masses. However, unlike all other 
methods, microlensing can in principle be used to measure the masses of objects without 
visible companions, in particular, single stars. At present, the Sun is the only single star 
whose mass has been directly measured with high precision. This was possible originally 
only because it has non-stellar, but nevertheless highly visible companions. 

In fact, the Sun is also the one single star whose mass has been accurately measured 
using gravitational lensing. While the original Eddington eclipse experiment was regarded 
at the time as a confirmation of Einstein's general relativity (Dyson, Eddington, & Davidson 
1920), general relativity is by this point so well established that this experiment can now 



1 The relationship between the observable parameters 6*e and f-E and the physical parameters M and 7r ro i 
is explained in Gould (2000). See especially his Fig. 1. 
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be regarded as a mass measurement of the Sun. Applying the same principle to 10 5 stars 
in the Hipparcos catalog, Froeschle, Mignard & Arenou (1997) were able to confirm general 
relativity (or alternatively measure the mass of the Sun) accurate to 0.3%. 

With the notable exception of the Sun, and despite the discovery of several thousand 
microlensing events as well as a decade of theoretical efforts to invent new ways to measure n E 
and 8 E , there have been just two mass measurements of single stars using microlensing. The 
problem is that while the microlens parallax 7r E has been measured for more than a dozen 
single lenses (Alcock et al. 1995, 2001; Mao 1999; Soszyhski et al. 2001; Bond et al. 2001; 
Mao et al. 2002; Smith, Mao k Wozniak 2002; Bennett et al. 2002; Smith, Mao & Paczyhski 
2003; Jiang et al. 2004), the angular Einstein radius 9 E has been measured for only five single 
lenses (Alcock et al. 1997, 2001; Smith, Mao & Wozniak 2003; Yoo et al. 2004; Jiang et al. 
2004). Though Alcock et al. (2001) and Jiang et al. (2004) each measured both 9 E and ir E for 
their events, respectively MACHO-LMC-5 and OGLE-2003-BLG-238, in neither case was the 
7Te measurement very accurate. Moreover, Gould (2004) showed that the microlens parallax 
measurement of MACHO-LMC-5 was subject to a discrete degeneracy. Nevertheless, Drake, 
Cook, & Keller (2004) resolved this degeneracy by a trigonometric measurement of 7r rc i. 
Gould et al. (2004) then combined the Drake et al. (2004) measurement of 7r re i and fi rcl with 
the original photometric data and additional high resolution photometry of the source to 
constrain the mass to within 17%. This is the most precise direct mass measurement of a 
single star (other than the Sun) to date. By comparison, the mass of the only other directly 
measured single star, OGLE-2003-BLG-238, is only accurate to a factor of a few. 

An et al. (2002) made the most precise microlens mass measurement to date, with 
an error of just 9%. However, the lens, EROS-BLG-2000-5, was a binary. In the future, 
the Space Interferometry Mission should routinely measure the masses of single stars both 
for stars in the bulge (Gould & Salim 1999) and for nearby stars passing more distant 
ones (Refsdal 1964; Paczyhski 1995; Salim & Gould 2000). Thus, at present, the direct 
measurement of single star masses (other than the Sun) remains a difficult undertaking. 

Here we present evidence that the microlensing event OGLE-2003-BLG-175/MOA-2003- 
BLG-45 is an excellent candidate for such a single-star measurement. This seems odd at first 
sight because, as we will show, 7i E is measured only to a factor of a few and 9 E is not measured 
at all. Hence it would appear difficult to derive any mass measurement, let alone a precise 
one. However, the event has the relatively unusual property that the lens itself is visible, 
and this makes a mass measurement possible. 

As discussed by Gould (2000) and in greater detail by Gould (2004), tt e is actually the 
magnitude of a vector quantity, ir E , whose direction is that of the lens-source relative motion. 
We first show that one component of tv e is extremely well determined, so that if its direction 
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could also be constrained, 7Te would also be well determined. Second, we show that the 
blended light for this event is almost certainly the lens. We outline how future space-based 
or possibly ground-based observations could measure /x rel , the vector lens-source relative 
proper motion (Han & Chang 2003). When combined with the very well determined £e for 
this event, this would yield #e through equation (1). At the same time, such a proper-motion 
measurement would give the direction of motion and so tightly constrain 7Te. 

2. Data 

The event [(RA,Dec) = (18:06:34.68, -26:01:16.2), (I, b) = (4.859, -2.550)] was initially 
discovered by the Optical Gravitational Lens Experiment (OGLE, Udalski et al. 1994) and 
was alerted to the community as OGLE-2003-BLG-175 through the OGLE-III Early Warn- 
ing System (EWS, Udalski 2003) on 2003 May 28. It was independently rediscovered by 
Microlensing Observations for Astrophysics (MOA, Bond et al. 2001) and designated MOA- 
2003-BLG-45 on 2003 July 6. It achieved peak magnification on HJD' =HJD-2450000 = 
2863.1 (2003 August 11). 

Observations were carried out by four groups from a total of eight observatories: OGLE 
from Chile, MOA from New Zealand, the Microlensing Followup Network (/zFUN; Yoo et 
al. 2004) from Chile and Israel, and the Probing Lensing Anomalies Network (PLANET, 
Albrow et al. 1998) from Chile, Perth, South Africa and Tasmania. OGLE made a total 178 
/ band observations from 2001 August 6 to 2003 November 10, of which 119 were during 
the 2003 season, using the 1.3m Warsaw telescope at the Las Campanas Observatory, Chile, 
which is operated by the Carnegie Institute of Washington. The exposures were 120 seconds 
and photometry was obtained using difference image analysis (Wozniak 2000). MOA made a 
total of 522 / band observations from 2000 April 12 to 2003 November 4, of which 303 were 
during the 2003 season, using the 0.6 m Boiler & Chivens telescope at Mt. John University 
Observatory in New Zealand. 

/zFUN monitoring of the event began on July 7. Observations were made at the 1.3m 
(ex-2MASS) telescope at Cerro Tololo InterAmerican Observatory in Chile using ANDICAM 
(DePoy et al. 2003) and at the Wise lm telescope at Mitzpe Ramon in Israel using the Wise 
TeK IK CCD camera. At CTIO, there were 210 observations in / band, from July 7 to 
October 29, and 11 observations in V, covering a similar period (July 9 - November 5). 
Exposures were 5 minutes each. Observations at Wise consisted of 12 in I band, covering 
the period July 8 to August 12, and 56 observations using a clear filter. The latter sampled 
the lightcurve densely just after peak, from August 12 to August 15. Photometry for all 
/xFUN observations was done using DoPHOT (Schechter, Mateo & Saha 1993). 
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PLANET observations of this event included: 52 observations in / band using the 0.9m 
telescope at CTIO, from August 11 to August 18; 80 observations in I band using the 0.6m 
telescope at Perth Observatory in Australia, from August 6 to November 2; 165 observations 
in R band using the Danish 1.54m telescope at La Silla, Chile, from June 4 to September 
1; 6 observations in / band using the South African Astronomical Observatory lm telescope 
at Sutherland, South Africa, on August 5; and 59 observations in / band using the Canopus 
Observatory lm telescope in Tasmania, from August 5 to September 21. The data reduction 
was done with the PLANET pipeline using PSF fitting photometry with DoPHOT. 

In fitting the lightcurve, we iteratively renormalized errors to obtain a x 2 P er degree 
of freedom of unity and eliminated points that were farther than 3a from the best fit. For 
the data sets (OGLE, MOA, /iFUN[Chile I, Chile V, Israel clear, Israel I], PLANET [Chile 
Danish R, Chile CTIO J, Perth, South Africa, Tasmania]), there were initially (178, 522, 
210, 11, 56, 12, 165, 52, 80, 6, 59) data points, of which (175, 515, 203, 11, 56, 12, 161, 51, 
76, 5, 51) were incorporated into the final fit, with corresponding renormalization factors 
(1.48, 1.179, 1.10, 1.00, 0.94, 0.81, 2.26, 1.06, 1.06, 2.69, 1.62). 

OGLE, MOA and /xFUN photometric data for this event are publicly available at 
http://bulge.astro.princeton.edu/~ogle/, http://www. roe. ac.uk/~iab/alert/alert. html and 
http://www.astronomy.ohio-state.edu/~microfun/. The data from all four collaborations 
are shown in Figure 1 together with a standard fit to the lightcurve, which shows strong 
residuals that are asymmetric about the peak. 



where F(t) is the observed flux, F s is the source flux, which is magnified by a factor A(t), and 
Fb is the flux from any stars blended with the source but not undergoing gravitational lensing. 
For point-source point-lens events, A(t) = A[u(t)], where u is the lens-source separation in 
units of 9 E and (Paczyhski 1986) 



3. Lightcurve fitting 



All microlensing events are fit to the functional form 



F(t) = F s A(t) + F h 



(4) 



A(u) 



u 2 + 2 



(5) 



u(u 2 + 4) 1/2 ' 



The event shows no significant signature of finite source effects, implying that equa- 
tions (4) and (5) are appropriate. However, it does show a highly significant asymmetry, 
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of the kind expected from parallax effects. We therefore fit for five geometric parameters 
(in addition to a pair of parameters, F s and F^, for each of the 11 observatory-filters com- 
binations). Three of these five are the standard microlensing parameters: the time of peak 
magnification, to, the Einstein crossing time, ts, and the impact parameter uo = u(t ). The 
remaining two are the micro lens parallax 7Te, a vector whose magnitude gives the projected 
Einstein radius, = AII/Ve, and whose direction gives the direction of lens-source relative 
motion. We work in the geocentric frame defined by Gould (2004) , so that the three standard 
microlensing parameters (t ,^E,w ) are nearly the same as for the no-parallax fit (in which 
7Te is fixed to be zero). 

The parallactic distortion of the lightcurve has a component that is asymmetric about 
the event peak, and one that is symmetric. The former allows a determination of 7Te,|| , 
the component of the parallax that is in the direction of the apparent acceleration of the 
Sun projected onto the plane of the sky at event peak. The symmetric distortion allows a 
determination of tte,±, the component perpendicular to 7Te,|| - The direction of 7Te,_l is chosen 
so that (7Te,||, 7Te,±) form a right-handed coordinate system. We fit for tte, however, as tte,n 
and 71e,e, the projections in the North and East directions (in the equatorial coordinate 
system), respectively. The error ellipse for these two 7Te parameters is highly elongated. To 
quantify this effect, we also calculate (7^,1,7^,2), the principal components of 7r E , as well as 
the position angle ip (north through east) of the minor axis of the error ellipse. The best-fit 
values thus obtained are shown in Table 1. However, microlensing lightcurve fits can suffer 
from several degeneracies. 

3.1. Degeneracies in the models 

Degeneracies arise when the source-lens-observer relative trajectory deviates from uni- 
form rectilinear motion but there is not enough information in the lightcurve to distinguish 
among multiple possible trajectories. We consider three types of degeneracies in our fits. 

3.1.1. Constant- acceleration degeneracy 

Since this is a relatively short-duration (£e ~ 63 days) event, the change in acceleration 
over this timescale is relatively small, and the fit is susceptible to the degeneracy derived by 
Smith et al. (2003) in the limit of constant acceleration. In the geocentric frame adopted 
in this paper, the additional solution is expected to have u' = —u , with the remaining 
parameters very similar to those of the original solution (Smith et al. 2003; Gould 2004). 
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That is, the lens passes on the opposite side of the source but otherwise the new trajectory 
is very similar to the old one. Table 1 shows that this is indeed the case. Moreover, the two 
solutions have almost identical x 2 - 

3.1.2. Jerk-parallax degeneracy 

Gould (2004) generalized the analysis of Smith et al. (2003) to include jerk and found an 
additional degeneracy whose parameters can be predicted analytically from the parameters 
of the original solution together with the known acceleration and jerk of the Earth at to- 
This prediction has been verified for both MACHO-LMC-5 (Gould 2004) and MOA 2003- 
BLG-37 (Park et al. 2004). We search for this potential alternate solution in two ways. First, 
we adopt a seed solution at the location predicted by Gould (2004) and search for a local 
minimum of the x 2 surface in the neighborhood of this seed. Second, we evaluate x 2 over 
a grid of points in the 7r E plane and search for any local minima. Neither search yields an 
additional solution. We note that for MACHO-LMC-5 (with timescale t E ~ 30 days) the 
two solutions have nearly identical x 2 , while for MOA 2003-BLG-37 (with t E ~ 42 days) the 
second minimum is disfavored at A% 2 ~ 7. It may well be that for events as long as OGLE- 
2003-BLG-175/MOA-2003-BLG-45 (t E ~ 63 days), the degeneracy is lifted altogether. 

3.1.3. Xallarap 

If the source is a component of a binary, its Keplerian motion will also generate accel- 
eration in the source-lens-observer trajectory. Like the Earth's motion, this is describable 
by the 7 parameters of a binary orbit. However, unlike the Earth's orbit, the binary-orbit 
parameters are not known a priori. Hence, while a parallax fit requires just two parame- 
ters, 7r E , (basically the size of the Einstein ring and the direction of the lens-source relative 
motion relative to the Earth's orbit), a full xallarap fit requires seven. This proliferation of 
free parameters may seem daunting but can actually be turned into an advantage in under- 
standing the event: if the full xallarap fit yields parameters that are inconsistent with the 
Earth's orbit, then this is proof that xallarap (rather than parallax alone) is at work. On 
the other hand, if the xallarap fit parameters are consistent with the Earth's orbit, this is 
evidence that parallax is the predominant acceleration effect. Of course, the latter inference 
depends on the size of the errors: if the xallarap parameters are tightly constrained and 
agree with the Earth's orbit, this would be powerful evidence. If the errors are very large, 
mere consistency by itself does not provide a strong argument. 
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We make two simplifications in our test for xallarap. First, instead of adding five 
parameters (to make the full seven), we consider a more restricted class of xallarap models 
with circular orbits. This eliminates two parameters, the eccentricity and the position angle 
of the apse vector. Hence only three additional parameters are required: the inclination, 
phase, and period. Second, rather than introduce additional free parameters into the fit, we 
conduct a grid search. 

We find that the data do not discriminate among the models very well. There is a large 
region of parameter space (including the Earth's parameters) that is consistent with the data 
at the 2a level. Only very short orbital periods, P < 0.2 yr are excluded. 

This exercise shows that, at least for this event, it is impossible to discriminate between 
parallax and xallarap from the lightcurve data alone. Hence, some other argument is needed 
to decide between these two possible interpretations of the acceleration that is detected in 
the lightcurve. 



4. Characteristics of the Blended Light 

We now argue that the blended light is most likely due to the lens. The key argument 
is astrometric: by measuring the centroid shift during the event, we show that the source 
and the blend are aligned to high precision and that the chance of such an alignment (if the 
blend were not associated with the event) is extremely small. In addition, the position of 
the blend on the color- magnitude diagram (CMD) shows it to be foreground disk star. 



4.1. Astrometry 

If we ignore the displacement of the positions of the images relative to that of the source 
(as justified below), then the position of the source-blend centroid of light, 6 C , is given by 
the flux- weighted average of the positions of the source and blend 2 : 

OUM , AF s [fl s + fx a (t - Q] + F h [0 h + ti h (t - Q] A*F s d s + F h h 

9e[m = afTTf, AF S + F b (6) 

where 6 S and #b are the positions of the source and blend at some fiducial time t*. n s and 
H h are the proper motions of the source and blend, and 0* is the centroid position at i*. For 



2 Alard, Mao, & Guibert (1995) calculate the shift in the centroid in the case of zero proper motion. Note 



that their equation (2) is in error and should read Ar c = r(l — /) 



1 Af+(\-f) 
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the time of maximum magnification (&* = to, A* = A max ), this equation can be rewritten 

e c [A(t)} = (0 b - e s )Z(A) + (A - l)n s W(A, t) + n c W(A, t) + O (7) 
where 6 is the centroid position at to, /x c = (fx s + rfj, b ), 

and 

We have introduced the parameter fi c instead of using \i h since the latter is highly 
correlated with \l s and the linear combination can be better constrained. For the OGLE 
data, r«2 (see Table 1). The quantities Z, W and A(t) are determined from the fit to the 
lightcurve. We fit the light centroid obtained from astrometry of 81 OGLE images taken 
both before and during the event to equation (7) and find, 

(0b - s ) North = 3.9 ± 7.6 mas, (0 b - # s ) E ast = -8.5 ± 10.5 mas. (10) 

The astrometric measurement errors of the individual points are assumed to be equal. Their 
amplitude is determined by forcing x 2 per degree of freedom to be unity. They are found to 
be 8 and 11 mas in the North and East directions, respectively. In Figure 2, we show this 
fit together with the data points plotted as A6 versus Z, where 

AO = e c - fj, s (A-i)w - n c w -e . (ii) 

Here the 6 C are the measured positions, while /x s , fi c , and 6 are the best fit parameters. 

Equation (10) shows that the source and blend have the same position within about 15 
mas. There are only 42 stars/arcmin 2 in this field that are as bright or brighter than the 
blend. The probability that one of them would lie within 15 mas of the source is therefore 
less than 10~ 5 , unless the star were related to the event. If the blend is related to the event, 
there are only three possibilities: (1) the blend is the lens, (2) the blend is a companion to 
the lens, or (3) the blend is a companion to the source. The last possibility is ruled out by 
the color-magnitude diagram (CMD), which shows that the blend lies in the foreground disk 
while the source lies either in or behind the bulge (see Fig. 3). While we cannot immediately 
rule out that the blend is a companion to the lens rather than the lens itself, we will show 
below that this hypothesis is ultimately testable. Moreover, even if the blend is a companion 
to the lens, most of the arguments of this paper remain unaltered. For the moment, we 
ignore this possibility and tentatively assume that the blend is the lens. 
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Before continuing, we note that neither of the proper-motion parameters is determined 
with high precision. We find 

/x s = (-28,-46) ± (64,89) mas yr" 1 and (/x s + rfi h ) = (5, 3) ± (5, 7) mas yr _1 . (12) 

This means, in particular, that the parameter of greatest physical interest (f-i h — jiij, which is 
a linear combination of these two fit parameters, can only be determined with a precision of 
about 100 mas yr -1 , far larger than its plausible value. Similarly, the astrometric errors are 
at least an order of magnitude too large to detect the motion of the image centroid relative 
to the source, which is why we ignore it in this treatment. 

Finally, we note that if the blend is either the lens itself or a companion to the lens, then 
one would expect the lens-source relative motion to be in the direction of Galactic rotation 
(roughly North by Northeast). This is because the blend lies in the foreground disk while 
the source lies in or behind the bulge. In fact, the parallax measurement shows that the 
lens-source relative motion is consistent with this direction (see Fig. 4). 



4.2. Mass and Distance of the Blend 



Independent of whether the blend is indeed the lens, we can obtain a rough estimate 
of the blend's mass and distance from its position on the CMD by making use of disk 
color-magnitude relation of Reid (1991), My = 2.89 + 3.37(V — I), together with the mass- 
luminosity relation of Cox (1999). This estimate necessarily involves a number of approxima- 
tions. First, the two relations just mentioned have scatter in them, which we ignore. Second, 
while the reddening could in principle be measured spectroscopically, no such measurement 
has been made. We therefore assume that the J-band extinction is related to the blend 
distance by 3 Aj = 0.5 (D^/kpc). Third, we must specify Ryj = Ay/E(V — /), the ratio of 
total to selective extinction. This is known to be anomalous toward the bulge, but while it 
varies somewhat from one bulge line of sight to another, the measured values lie consistently 
near R VI ~ 2.1 (Popowski 2000; Udalski 2003; Sumi 2004). We therefore adopt this value. 

Fourth, we must estimate the apparent color and magnitude of the blend. While in 
principle the most straightforward step, under present circumstances this is actually the 
most uncertain. The flux of the blend is a parameter of the fit to the lightcurve. To 



3 Extinction is proportional to exp 



JL 



where z is the height above the Galactic plane, R is the 
distance, along the plane, from the Sun to the blend, and Zh and R s are the dust scale height and scale 
length, respectively. R has its origin at the Sun and increases towards the Galactic Center. Since z = Rsinb, 
for Zh ~ 130 pc, R s ~ 3 kpc and b = 2.5° the two terms inside the exponent almost cancel. 
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determine a color, we must have two such fluxes and so use the /zFUN Chile photometry 
since this is the only one of our observatories with data in two photometric bands. However, 
it is known that the /xFUN Chile photometry contains additional blended light relative to 
the OGLE photometry. First, the ratio of fit parameters, r = F^/F S) is greater for /iFUN 
Chile I (r ~ 2.60) than for OGLE / (r ~ 2.01) despite the fact that the passbands are very 
similar. Second, OGLE photometry identifies additional sources in the neighborhood of the 
source that /iFUN photometry does not identify and that therefore must be included in the 
/iFUN blend. If the colors of these extra sources were known, they could just be removed 
to find the color as well as the magnitude of the OGLE blend. Unfortunately they are not 
known. For the purposes of this estimate we assume the color of the blend is also the color 
of the lens. However, the fact that the better-determined OGLE ratio r ~ 2.01 is less than 
r ~ 2.60 for /iFUN indicates that the lens may be fainter by 0.28 mag than its /iFUN I 
magnitude. We therefore use this corrected value for our estimate. Finally, the CMD has 
not been directly calibrated to standard bands. The OGLE I fluxes are calibrated to within 
a few tenths, and by identifying the OGLE and /iFUN J, we can therefore approximately 
calibrate the ordinate of the /iFUN CMD. We then determine Aj = 2.0 of the clump from 
the calibrated / of the clump and the known dereddened magnitude of the clump, I Q = 14.32 
(Yoo et al. 2004). We then estimate E(V — I) = 1.82 using R V i = 2.1 and so from the known 
dereddened color of the clump, (V — I)o = 1.00, calibrate the abscissa. Clearly, the very 
complexity of this approach as well as the sheer number of approximations leaves something 
to be desired. Nevertheless, as we are interested only in rather crude mass and distance 
estimates for the blend, it will suffice. We estimate, 

J b = 16.48, (y-/) b = 1.67. (13) 

To carry out our calculation, we consider trial stars as a function of blend mass, M h . 
For each mass, we obtain an absolute magnitude My using Cox (1999) and then a color 
{V — /)o,b using the Reid (1991) color-magnitude relation. This gives a selective extinction 
E(V - I) = (V - I) h - (V - J) ,b and so an extinction A T = E(V - I){R V i - 1), and 
hence a distance = Aj/(0.5 mag/kpc). From this, we obtain a predicted / magnitude, 
/ = Mi + 5 log-Db/10pc + Aj. These predictions are shown in Figure 5 where they are 
compared to the observed magnitude lb = 16.48. From this comparison, we obtain 



M b = O.75M , D b = l.lkpc. 



(14) 
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4.3. Microlens Parallax and Proper-Motion Predictions 

If we identify the blend with the lens, and assume that the source lies at D s ~ 10 kpc, 
we obtain 7r rc i = n\ — it s = 0.91 mas — 0.10 mas = 0.81 mas. Assuming M\ = Mb = O.75M , 
and substituting these values into equations (2) and (3) yields, 

tt e = 0.36, 9 E = 2.2 mas. (15) 

We may now ask if these values, which are derived from the photometrically-determined 
characteristics of the blend, are consistent with what is known about the microlensing event. 
From Figure 4, we see that the predicted parallax, 7Te = 0.36, is consistent with the value 
observed at the ~ 1 a level. Combining the Einstein radius 6>e = 2.2 mas with the event's 
measured timescale t& = 63 days yields a proper motion 

/z re i = 7^ = 13 masyr -1 . (16) 

The only hard information we have on /i rcl comes from the lack of finite source effects, 
which puts a weak lower limit on the Einstein radius, 9 E > 9^/uq. Here 9* is the source radius 
and Mo — 0.05 is the impact parameter. Using the standard method to infer the angular 
source size from the instrumental CMD (Yoo et al. 2004), we find 9* = 3.8/ias. Hence 9 E > 
76/ias, and /i rc i > 9 E /t E — 0.5 masyr -1 . Equation (16) easily satisfies this limit. However, 
the proper motion in equation (16) is somewhat higher than the typical (/i rc i ~ 7 masyr -1 ) 
proper motion that would be expected for a disk lens moving with same rotation velocity 
as the Sun and seen projected against a star with some random motion in the bulge. But, 
given that the blend is so close (D h ~ 1.1 kpc), the peculiar motions of the Sun and the 
blend relative to the mean disk rotation may both contribute significantly to fj, TC \. Finally, 
the measurement of the blended motion /x c = fj, s +r/j, b = (5, 3)±(5, 7) masyr -1 (see eq. [12]), 
also places indirect constraints on /x rel = [fi c — (r+l)/it s ]/r. That is, since r ~ 2, equation (16) 
implies |3/x s — fi c \ = 24 masyr -1 . This constraint is not easily satisfied unless either fi s is 
anomalously fast or fi s is anti-aligned with /x c (and so /j, b ). That is, /x s • /x c < 0. However, 
the latter option is quite plausible. The source could be retrogressing, as would occur if it 
were in the far disk, and as would be consistent with its position somewhat below the clump 
in the CMD. In that case, the relative proper motion could be high without requiring rapid 
motion of the centroid of light. Hence, the proper motion obtained in equation (16) is not 
unreasonable. 

In brief, all the available evidence is consistent with the hypothesis that the blend is the 
lens or a companion to the lens. In either case, this opens the possibility that the lens mass 
can be precisely determined by measuring the proper motion of the blend. There would still 
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remain the question of whether the mass that was measured was that of the blend, in which 
case this measurement could be compared with more accurate photometric and spectroscopic 
measurements than have been obtained to date. We return to this question in § 6. 



4.4. Uncertainty in Lens Mass Estimates 

The mass M of the lens is given by 



M=-^ = ^. (17) 
KTTe kite 



The uncertainty in the mass estimate is therefore 



2 / \ 2 / \ 2 



In microlensing events in general, the fit to the lightcurve typically produces a tight con- 
straint on t E , so <7t E /tE is usually small. Microlensing parallax is usually less well-determined. 
In this event, for example, the fit constrains very well the component of 7Te that is parallel to 
the direction of acceleration at the event peak but the other component is only very poorly 
constrained. We refer to these components 4 as 7Te n and 7Te,±. This situation is depicted in 
Figure 6, where the solid line I shows the direction of the long axis of the error ellipse for 7Te 
(see Fig. 4 for comparison), and the dashed lines parallel to I and at a distance a WE ., indicate 
the uncertainty. A hypothetical measurement of the direction of the relative proper motion 
is shown in Figure 6 as the line to, with its associated uncertainty indicated by the dotted 
lines on either side. The magnitude of 7Te is then given by the distance from the origin to 
the point of intersection of lines / and to (OA in Fig. 6). 

The uncertainty in 7r E thus has two contributions, one from the width a nE |( and the 
other from a fl± . The errors are not correlated and so may be added in quadrature. The 
fractional uncertainty in 7r E is therefore given by 

fey-sH^ftH' <i9) 



4 For short events, < yr/27r, the short axis of the error ellipse should line up with the direction 
of acceleration (Gould, Miralda-Escude, & Bahcall 1994) and this prediction has been confirmed to high 
precision for two short events (Park et al. 2004; Jiang et al. 2004). For OGLE-2003-BLG-175/MOA-2003- 
BLG-45 the acceleration position angle is 87.1°. So while the minor axis of the error ellipse for the uo > fit 
is aligned with the direction of acceleration, that for the uo < fit differs by 2.4°. We ignore this difference 
in this section, and refer to the principal components as "tte.h" and "7Te,_l". 
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where 7 = 90° — (a + (5) and the angles a, (3 and 7 are as defined in Figure 6. 

For this event, the fit to the lightcurve gives (a tE /t E ) = 0.012 and (av B n/vr E ,j|) = 0.056, 
both small. The error in the mass will therefore be dominated by the uncertainty in the 
proper motion measurement, unless it is very accurate. If (<7 Mrel //i re i) = (cr M± / '// re i) and 
(a + P) w 45°, a 10% determination of // re i would give a 16% measurement of the mass. 



5. Measuring the Proper Motion 

The actual determination of the lens mass depends on an accurate measurement of the 
lens-source relative vector proper motion: the magnitude is required to determine 6>e = /irei^E, 
while the direction is required to determine 7Te (see § 4.4). 

There are two ways in which /x rel may be observationally determined. It will be possible, 
by waiting long enough, to resolve the source and the blend into separate objects in ground- 
based images. However, at the estimated rate of 13masyr _1 that might take of the order of 
a decade. 

The other method uses the higher resolution of space-based imaging and the significantly 
different reddening of the source and the blend. The blend is on the reddening sequence on 
the CMD, which indicates it lies in front of most of the dust column to the source. The 
source, which lies in or behind the bulge, is therefore much more reddened than the blend. 
Hence, by imaging in a blue band such as U, in which the source is expected to suffer high 
extinction, and in a red band such as J, in which it suffers much less extinction, it would be 
possible to measure the separation between the source and the blend from the offset of their 
centroids in the two bands even before they are separately resolved. The true separation AO 
is given in terms of the separation AOm in the U and I band centroids by 

1 1 



M 1lT»-ITmJ Ae "' (20) 

where {F a /F b )u and (F s /F b )j are respectively the source/lens flux ratios in U and I. The 
proper motion is then simply /x rel = AO /At, where At is the time interval between the 
event peak and the epoch of the observation. With the resolution of the Hubble Space 
Telescope, such a measurement should be possible 3 years after the event peak. Note that 
while (F s /Fb)i is known from the microlensing event itself, the determination of (F s /Fb)u 
requires a bit more work. First, (F s + Fb)u is measured directly in the followup observations. 
To find F Sj u, one should first note that (V — I) s is well determined from the microlensing 
fit. Hence, (U — I) s will be very similar to the (U — I) of other clump stars with the same 
(V — I) as the source. One can evaluate the error in this determination from its scatter 
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when applied to other clump stars. This procedure yields F St u, and so (together with the 
total-flux measurement) (F s /F b )jj. 

Applying this method to blue and red plates from the 1950 Palomar Sky Survey, we 
find a shift A0 BR = (0.30 ± 0.17,-0.23 ± 0.17) arcseconds, with the red centroid to the 
north-east of the blue. While this direction agrees with the relative proper motion expected 
for a disk lens, the magnitude is consistent with zero and provides no useful constraint. 

5.1. Geocentric versus heliocentric frames 

Though we have been working in the geocentric frame in this paper, the lens-source 
relative proper motion fj, iel is measured in the heliocentric frame. Hence, to measure the 
mass M = /i rc i^E,hci/ K7r E) the Einstein timescale t E = tE, g eo that we obtain from fitting 
the lightcurve must be transformed to the heliocentric frame, in which it is not as well 
determined. The geocentric and heliocentric timescales are related by r E = tE,geo^geo — 
^E.hei^heb where r E is the projected Einstein radius and v is the projected velocity in the 
appropriate frame. The transformation of the projected velocity to the heliocentric frame is 
accomplished using the known geocentric velocity of the Sun at event peak: Vh e i = v geo — v . 
Since v oc 1/(7Te£e) and t E is well-determined in the geocentric frame, the uncertainty in v geo 
is dominated by the uncertainty in n E shown by the la contours in Figure 4. As shown 
in Figure 7, the corresponding uncertainty in the ratio iE,hei/*E,geo = v gco /vhc\ is about 3%, 
almost independent of the actual direction of the proper motion. 

6. Distinguishing the Lens and Blend Hypotheses 

If the mass of the lens is measured from the proper motion of the blend, it will still not 
automatically be known whether the lens is the blend or is a companion to it. Here we show 
that the proper-motion measurement itself can help distinguish these hypotheses. 

There is no sign of binarity in the well-sampled lightcurve. This provides a lower limit 
to the binary separation if the lens system is a wide binary, or an upper limit if the system 
is a close binary, as follows. 

At large separations, the companion would induce a caustic in the magnification profile 
of full width £ = 4g~ 1 rf -2 , where q is the lens/blend mass ratio, d6 E is their separation, and 
£9 E is the angular width of the caustic. Since the source clearly did not traverse a caustic, 
£ < 2v / 2mo ) and indeed, detailed fitting would provide somewhat tighter constraints (see Fig. 
1 of Gaudi & Gould 1997). Hence, d > (quo/^/2)- 1 / 2 = 5.3<T 1/2 . From equation (15), the 
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Einstein radius associated with the blend is expected to be 6e = 2.2 mas. The lens Einstein 
radius would be smaller by q x l 2 . Hence, the separation of the lens from the putatively distinct 
luminous blend must have been at least 

|0j _ Q h \ = dq 1/2 9 E:h > 12 mas (21) 

at the peak of the event. 

This separation is already of order the measurement errors from the OGLE astrometry 
(see eq. [10]). If future proper- motion measurements are taken at multiple epochs, they 
should be able to determine whether the blend-source relative motion points back to a 
common position at the time of the event, or whether the two were separated by at least the 
lower limit from equation (21). Such a measurement would therefore be able to determine 
whether the blend was the lens, or a companion to the lens. 

If the system is a close binary (either one of the components is unseen, or both are 
visible but unresolved), then d < (q 1 ^ 2 + q~ l ^ 2 ){uQ/ ' v /2) 1 / 2 . For small q this limit can be 
approximated as d < 1/(5.3 q 1 ^ 2 ). Since v 2 = M/dr-E, where v is the orbital velocity of the 
secondary, that implies v > 1.6g 1 / 4 (M/M ) 1 / 2 f e , where we have used 9e = 2.2 mas and 
ix\ ~ 1 kpc. The radial velocity of the primary is then qv > 42g 5//4 km s^ 1 . This should be 
detectable unless the companion is substellar (g < 0.1). 
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Table 1. Best-Fit Parameters 





Uo 


> fit 




< fit 


Parameter 


Value 


Uncertainty 


Value 


Uncertainty 


t (days) 


2863.1116 


0.0065 


2863.1119 


0.0059 




0.0546 


0.0008 


-0.0547 


0.0008 


t E (days) 


62.7894 


0.8874 


63.2196 


1.1297 


7I"E,N 


0.1108 


0.7296 


0.2124 


0.3463 


7TE,E 


n 1 

U.lOUo 


U.UooO 


U.14oo 


U.UooO 


7I"E,1 


0.1658 


0.0090 


0.1674 


0.0090 


7TE,2 


-0.1024 


0.7306 


-0.1977 


0.3478 




87.1° 




84.7° 






2.0057 


0.0020 


2.0042 


0.0018 


(F fc /F s ) 2 


2.5965 


0.0027 


2.5950 


0.0024 


(F b /F s ) 3 


2.6459 


0.0029 


2.6439 


0.0026 


(F b /F a ) 4 


3.4187 


0.0071 


3.4172 


0.0067 


{F b /F s ) 5 


2.7371 


0.0109 


2.7335 


0.0106 


(F b /F s ) 6 


2.8999 


0.0048 


2.8974 


0.0045 


(F b /F s ) 7 


3.9967 


2.7194 


3.9964 


2.7173 


(F b /F s ) 8 


2.6091 


0.0030 


2.6069 


0.0028 


(F b /F s ) 9 


7.0954 


0.0317 


7.0921 


0.0301 


(F b /F s ) w 


4.2482 


0.1047 


4.2433 


0.1039 


(F b /F.) u 

x 2 


3.4974 


0.0050 


3.4946 


0.0045 


1296.6528 




1295.9659 





Note. - Observatory/filter combinations for the ratios (F b /F s ): 
l=OGLE 7, 2 =y uFUN Chile 7, 3=MOA 7, 4= y uFUN Wise 7, 5=PLANET 
CTIO 7, 6=PLANET Perth 7, 7=PLANET SAAO 7, 8=PLANET Tas- 
mania 7, 9=mFUN Chile V, 10=^FUN Wise clear, 11=PLANET Danish 
R. 
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Fig. 1. — Fit to the lightcurve in the nonparallax (left panel) and the parallax (right panel) 
models. Data from different collaborations are color-coded as follows: //FUN, red; MOA, 
green; OGLE, black; PLANET, blue. Symbols indicate filter bands: J, open circles; V, 
inverted triangles; clear, filled circles; R, upright triangles. Data have been binned by 
night. The residuals from the nonparallax fit show an asymmetry about the event peak; the 
asymmetry disappears when the lightcurve is fit for parallax. The inset shows a zoomed-in 
view of the peak. 
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Fig. 2. — Measurement of the separation of the microlensed source and the unlensed blended 
light from OGLE astrometry. As described in § 4.1, the blend-source centroid should be 
at O c [A(t)] = (0 b - e s )Z{A) + (Afx s + rn h )W(A,t) + where r = F h /F s , W(A,t) = 
(t - t ) /(A + r), and Z(A) = (A mSLX - A)r/[(A max + r)(A + r)]. The two panels show the 
residuals, using the best-fit values of fi s , fi h and ^o, in the East and North directions but 
with the (#b — ^s)Z(A) term removed. Hence the slopes of the linear fits to these residuals, 
($b — $s)North = 3.9 ± 7.6 mas and (#b — ^s)East — —8.5 ± 10.5 mas give estimates for the 
source-blend separation. This separation is consistent with zero within small errors. 
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Fig. 3. — Approximately calibrated color-magnitude diagram of stars within a 6' square of 
OGLE-2003-BLG-175/MOA-2003-BLG-45. See § 4.2 for details on the calibration procedure. 
The source (filled circle) lies a magnitude below the center of the bulge clump (triangle) and 
therefore lies within or behind the bulge. The open circle shows the position of the blend 
based on /xFUN photometry, and the cross shows the position after the correction based on 
OGLE photometry. The blended light lies along the "reddening sequence" of foreground disk 
main-sequence and turnoff stars, which suffer less extinction than the bulge and therefore 
appear relatively bright and blue in the diagram. Based on this diagram and arguments 
given in § 4.2, we conclude that the blended light comes from a star of mass M h = 0.75 M 
and distance D h = 1.1 kpc. 
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Fig. 4. — Left panel: la contours of the x 2 surfaces, projected onto the (7Te,e, ^e,n) plane, 
for the u > and uq < fits. The arrow shows the direction of Galactic rotation. The 
axis ratio of the Uq > "ellipse" is 26 : 1 and that of the Uq < one 24 : 1. Though only 
one component of 7Te is well-constrained, this figure shows that the direction of lens-source 
relative proper motion is consistent with the direction of Galactic rotation, which supports 
the hypothesis that the lens is a foreground disk star. The dotted arcs are parts of a circle 
of radius 7r E = 0.36, the predicted microlens parallax. The inset shows the contours with 
the horizontal scale expanded. Right panel: The solid contours are the same ones shown 
in the left panel, and are in the geocentric frame. The contours in dashed lines show them 
transformed to the heliocentric frame. 
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Fig. 5. — Expected apparent / magnitude vs. mass, as explained in § 4.2. The observed 
I = 16.48 magnitude of the blend is marked with a dotted line. The intersection of the two 
curves gives our estimate of the mass of the blended light source, M b = 0.75 M Q . 



o 



7tE,East 

Fig. 6. — Even when only one component of the parallax 7Te is well-determined, it may 
be possible to constrain the magnitude he if the source-lens relative proper motion /z rel 
is measured. In this figure it is assumed that the (7Te,||, tte,±) axes are rotated from the 
(7TE,East> 7r E ,North) axes by an angle a, that 7Te,|| is well-constrained and that 7Te,± is not con- 
strained at all. The solid line I is the (very) long axis of the error "ellipse" , which in this case 
is a strip of width 2PP' = 2a 7TE .. . The solid line m shows the direction of relative proper 
motion, and the dotted lines on either side indicate the uncertainty in the direction. The 
magnitude of 7r E corresponds to the length of OA. The uncertainty in the length of OA 
results from the uncertainty PP' in 7Te,|| and the uncertainty AA' in the direction of /x rel . 
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Fig. 7. — Ratio of Einstein timescales in the heliocentric and geocentric frames, evaluated 
along the contours shown in Figure 4 for the Mo > (solid line) and the uq < (dashed line) 
fits. For a given direction of lens-source relative proper motion (3 = tan _1 (7r E Noi . th /7r E)East ), 
(/3 can be inferred from the direction of /u rel ) the width of the la contours corresponds to a 
variation in tE,hci/^E, g co of about 3%. 



